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ABSTRACT 


We investigate the utility and robustness of a new statistic, ujj (rg), for analyzing Baryon 


Acoustic Oscillations (BAO). We apply oOf (vc), introduced in Xu et al. (2010), to mocks and 
data from the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey 
(BOSS) included in the SDSS Data Release Eleven (DRll). We fit the anisotropic clustering 
using the monopole and quadrupole of the uif (vc) statistic in a manner similar to conventional 
multipole fitting methods using the correlation function as detailed in ( Xu et al.|2012| ). To test 
the performance of the uif (vc) statistic we compare our results to those obtained using the 
multipoles. The results are in agreement. We also conduct a brief investigation into some of 
the possible advantages of using the uif (vc) statistic for BAO analysis. The uif (vc) analysis 
matches the stability of the multipoles analysis in response to artificially introduced distortions 
in the data, without using extra nuisance parameters to improve the fit. When applied to data 
with systematics, the uif (vc) statistic again matches the performance of fitting the multipoles 
without using nuisance parameters. In all the analyzed circumstances, we find that fitting the 
uif (vc) statistic removes the requirement for extra nuisance parameters. 

Key words: cosmology:large-scale structure of universe 


1 INTRODUCTION 

Baryon Acoustic Oscillations (BAO) are a powerful tool to break 
degeneracies in measurements that rely solely on the cosmic mi¬ 
crowave background (CMB). For example, CMB measurements 
alone cannot obtain satisfactory constraints on the Hubble constant, 
Ho, and the matter density, Qrn ( [Hinshaw et al.|201^ . 

BAO arise from acoustic oscillations before recombination; 
their effects are imprinted on the galaxy distribution in the universe 
( [Eisenstein & Hu|1998| ). The scale of these oscillations can be used 
as a standard ruler to provide distance measurements to various red- 
shifts. These measurements are sensitive to the value of cosmolog¬ 
ical parameters and provide independent observations which break 
CMB degeneracies. 

With the most recent galaxy surveys, it is possible to measure 
the BAO feature in both the line-of-sight and transverse directions. 
This constrains both the Hubble parameter, H (z), and the angular 
diameter distance. Da {z), at redshift z. Because the BAO feature 
is sensitive to cosmology, the feature shifts in the line-of-sight and 
transverse directions based on the fiducial cosmology used during 
analysis. By applying the Alcock-Paczynski test, we can use these 
measurements to constrain cosmology jAlcock & Paczynski|1979| ). 

Over the past few years, several statistics have been proposed 
for analyzing the BAO feature in galaxy surveys. In particular, the 
multipoles fitting method has been applied to data from the SDSS- 
III BOSS Data Releases to produce consistent constraints on H (z) 


and Da (z) ( |Xu et al.|2012| ). In this paper, we focus on a statistic, 
(^c), which may provide several advantages over ones in current 
use (Xu et al. 2010| ). We apply the Ui (vc) statistic to data and mock 
catalogs from the BOSS Data Release 11 and provide comparisons 
to constraints produced by the multipoles method. We also perform 
a brief investigation into the possible advantages of working with 
the U£{rc) statistic. 

The remainder of this paper is structured as follows. In Sec¬ 
tion!^ we describe the data and mocks used in this paper. In Sec¬ 
tion we give a brief introduction of the theory behind the fitting 
method used. In Section|^we present the details of the fitting model 
used in our analysis. In Section|^we present our results from apply¬ 
ing the Ui (tc) statistic to mock catalogs and data from the SDSS- 
III Data Release 11 (DRll). In Section]^ we describe our investi¬ 
gation into the advantages of using the Ui (fc) statistic. Finally, we 
conclude in Section |7] 


2 DATA AND MOCKS 
2.1 Data 

We use data included in Data Release 11 (DRl 1) of the Sloan Dig¬ 
ital Sky S urvey (SDSS; [York et al.| ( [T000| )). Using the 2.5m Sloan 
Telescope ( |Gunn et al.|2006| ) at Apache Point Observatory in New 
Mexico, SDSS I, II (|Abazajian et al.|20091), and III (|Eisenstein et 
|^[20TT] ) used a drift-scanning mosaic CCD camera ( [Gunn et al"] 
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|1998| ) to image 14,555 square degrees of the sky in five photomet¬ 
ric bandpasses ( [Fukugita et al.|(T9^ ; [Smith et al.| ( [200^ ; |Doi et| 
[n?] ( |2010| )) to a limiting magnitude of 22.5. The data was passed 
through pipelines designed to perform astrometric calibration ( |Pier| 
|et al.|2003] l, photom etric reduction {[Lupton et al.|2001| l, and pho¬ 
tometric calibration ( [Padmanabhan et al.|2008| l. All of the imaging 
data was re-processed as part of SDSS Data Release 8 (DR8;|A^ 
|hara et alTl ( |2011| ). 

The Baryon Oscillation Spectroscopic Survey (BOSS) itself 
is designed to obtain spectra and redshifts for 1.35 million galax¬ 
ies over a footprint of approximately 10,000 square degrees with a 
redshift completeness of over 97 percent. The targets are selected 
from SDSS DR8 imaging and assigned to tiles of diameter 3° us¬ 
ing a target density adaptive al gorithm ([Blanton et al.|200^ . Red- 
shifts are derived from spectra ( [Bolton et al.|2012|) obtained using 
the doubled-armed BOSS spectrographs ( Smee et al.|2013j t. For a 
summery of the survey, one may consult Eisenstein et al.|f2011j ). 
[Dawson et al.| ( |201^ also provides a full description. 


2.2 Simulations 

We use 600 SDSS III-BOSS DRll PTHalos mock galaxy catalogs 
with the same angular and radial masks as the survey data to com¬ 
pute the sample covariance matrices used in our analysis. These 
mock catalogs, provided by [Manera et al.[p013[ ), are generated at 
^ = 0.55 in boxes of side length 2400/i“ Mpc with 1280^ dark 
matter particles. 

When generating the mocks, [Manera et al.[ ( [2013[ ) first used 
second-order Lagrangian perturbation theory (2LPT) to create a 
matter density field. Then, the halos were identified using a friends 
of friends halo finder. Masses and linking lengths were calibrated 
to N-body simulations, and the halos were populated using a Halo 
Occupation Distribution calibrated to match the observed cluster¬ 
ing in CMASS catalog ( [White et al.|2011] l. 


that parameterizes the BAO shifts. The parameters can then be mea¬ 
sured by fitting the model to our data. As in [Xu et ^ ( [2012[ ), we 
parameterize the isotropic shift using a and the anisotropic signal 
using e. These parameters are defined by: 

Dv{£)lrs 


1 + e — 


Hfid ( 2 :) Da, fid ( 2 ^) 
H{z) Da (z) 


1/3 


( 2 ) 


Here, Ts denotes the sound horizon or the BAO scale. The subscript 
fid denotes the fiducial cosmology (see Section [4TJ which we as¬ 
sume when making our measurements. Dy {z) is the spherically 
averaged distance to redshift z. It is defined by: 


Dv (z) = 


{l + z)^D\ (z) 


cz 


1/3 


H{z) 


( 3 ) 


Note that if there is no isotropic shift, then a = 1 and if there 
is no anisotropic warping, e = 0. 

Alternatively, we can also parameterize using the shift parallel 
to the line-of-sight, a\\, and the shift perpendicular to the line-of- 
sight , a±: 


Da {^Z^ Ts,fid 

Daju {z) Ts 


(4) 


a± = 


Hfid ( 5 :) Tg^fid 


(5) 


H {z) Tg 

The relations between the two parameterizations are given by: 


2/3 1/3 
— a / af 


( 6 ) 


1 + e — 



(7) 


3 BACKGROUND AND THEORY 
3.1 Background and Parameterization 

There are two main effects that contribute to anisotopic cluster¬ 
ing: redshift-space distortions and the adoption of an incorrect cos¬ 
mology when calculating galaxy separations. Redshift-space dis¬ 
tortions are generated by the peculiar motions of galaxies within 
clusters (Finger-of-God effect) and by larger scale fiows of galax¬ 
ies into overdense regions ( [Kaiser These effects distort our 

measurements of line of sight separation between galaxies and af¬ 
fect the shape of the correlation function. 

Adopting an incorrect cosmology distorts our measurements 
of galaxy separations in both the line-of-sight and transverse di¬ 
rections. Line of sight separations are computed using the Hubble 
parameter H (z) while transverse separations are computed using 
the angular diameter distance Da (z). However, both H (z) and 
Da {z) are dependent on cosmology. As a result, adopting an in¬ 
correct cosmology gives a different clustering signal and BAO sig¬ 
nal in the line-of-sight direction versus the transverse direction ( [Xu[ 
[et al.|201^ . In addition to these anisotropic effects, the BAO fea¬ 
ture can also be shifted isotropically as a result of assuming an in¬ 
correct cosmology. By applying the Alcock-Paczynski test, we can 
use these shifts to constrain our cosmology ([Alcock & Paczynski[ 
[T?7^ . 

To measure the anisotropy we construct a clustering model 


3.2 Clustering Estimators 


Most clustering estimators for measuring the BAO feature require 
either the computation of the 2D correlation function, ^ (s, /x), or 
the power spectrum, P (k) ( [Anderson et al.[[2014b|). The two es¬ 
timators we investigate here. M ultipoles ([Xu et al.|2012[ ) and the 
Ui (tc) statistic Ee eF^(20T0t ; rBlake et 31.^201!^ are based on 
analysis of the 2D correlation function. However, working with the 
full 2D correlation function is impractical as it requires a far larger 
quantity of mock catalogs than are currently available. To avoid this 
problem, we perform our analysis by compressing the correlation 
function into a small number of angular moments. 


3.2.1 Multipoles 

Since the primary focus of this paper is the ujg (vc) statistic, we will 
only present a summary of the multipole analysis. The interested 
reader may consult [Xu et al.[ ( |2012[ ) for a more detailed discussion. 
First we define our coordinates: 

22,2 

r =r||+r^ (8) 

jji — cosO — — (9) 

r 

Here, r is the separation between two galaxies, r|| is the separa¬ 
tion in the line-of-sight direction, and ri_ is the separation in the 
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transverse direction. 0 is the angle between a galaxy pair and the 
line-of-sight. 

Then, a and e are defined by: 

r\\=a{l + ef r\\ ( 10 ) 

r'±= a{\ + e)“^ rx (11) 

where the primed coordinates denote the true cosmology space and 
the unprimed the fiducial cosmology space. 

The multipole analysis measures the Legendre moments of the 
2D correlation function: 

6 W = “^^2 ^ ^ *^12) 

where Lg (/i) is the ^th Legendre polynomial. 

In the multipole analysis, we focus on the monopole and 
quadrupole since the infiuence of higher order terms is negligible. 
Note that the odd order multipoles are zero due to symmetry. 


3.2.2 The ug {tc) Statistic 

As in Xu et al.| ( |2010| ), we define uji (vc) as the redshift space cor¬ 
relation function, {r, /i), convolved with a compact and compen¬ 
sated filter Wi (r, /x, rc) as a function of characteristic scale Fc'. 


!^e {rc) j (^> A*) (^> A*’ '^c) 

= j (r, jj) Wi (r, Tc) Lt {n) (13) 

= j dr (r) Wi (r, rc) 


where we have taken advantage of the orthogonality of the Legen¬ 
dre polynomials and set Wi (r, /x, rc) = Wg (r, rc) Lg (/x). 

Following [Padmanabhan et al^ ( |2007| ) and |Xu et "^ ( 120101 ) we 
define a smooth, low order, compensated filter independently of t. 


W (x) = {2xf {1 - xf ^ (14) 

where x = (r/rc)^. 

This choice of filter gives the Ui (rc) statistic several advan¬ 
tages. By design, Ui (rc) probes a narrow range of scales near the 
BAO feature and is not sensitive to large scale fiuctuations or to 
poorly measured or modeled large scale modes ( |Xu et al.]|2010| ). 
Further, Ui (rc) is blind to the constant ^ monomial, which is the 
integral constraint term that Ui (rc) is designed to avoid. Any other 
monomials in ^ also produce the same monomial in Ui (rc). 


3.3 Modeling the correlation function 

A model for the correlation function is required in order to per¬ 
form our fitting analysis with our two estimators. We start with a 
nonlinear power spectrum template: 

P(fc,At) = (l + /3/i")V(fc,At,E,)Pd„(fc,/i,Ex,S||) (15) 

where the (l + y5/x^)^ term corresponds to the Kaiser model 
( |Kaiser| |1987| ) for large scale redshift space distortions. The 
F {k, /X, Ss) term is the streaming model for the Finger-of-God ef¬ 
fect ( [Peacock & Dodds |1994| ): 


where Eg is the streaming scale. Pdw (^,/x, E||) is the “De- 
Wiggled” power spectrum adopted by jXu et ^ ( |2012j ) and jAnder^ 
[son et al.] ( |2014a| ), among others, as a template for the non-linear 
power spectrum. E^ and E|| serve to parameterize the Gaussian 
damping of BAO due to nonlinear structure growth (jEisenstein etj 
|al.|2007| l. 

We produce the multipoles of the correlation function by first 
decomposing the 2D power spectrum given in Eq.fTS] 

Pi {k) = y, ^ 

We then transform to configuration space: 

(r) = J (fc) je {kr) (18) 

where {kr) is the ^th spherical Bessel function. 


3.4 Covariance Matrix 


In this work, we apply corrections to the covariance matrix as sug¬ 
gested in jPercival et al.|p013| ). Further, since we estimate the co- 
variance matrix from mock catalogs which parallel BOSS observa¬ 
tions by splitting the observations into the Northern Galactic Cap 
(NGC) and the Southern Galactic Cap (SGC), we must adjust the 
standard computation for the covariances accordingly. For a de¬ 
tailed description of our computation, the reader may consult [Wr^ 
[gas Magana et al.| ( [2013| ). 

The ij\h entry in the covariance matrix, C, used in our fitting 
analysis is computed using the 600 PTHalos mocks: 

-^300 

m>300 


where ^ = Ko, ^ 2 ], is the average over all the mocks, and is 
the correlation function for the mth mock. The first 300 mocks cor¬ 
respond to the NGC and the rest correspond to the SGC. When 
computing C for the ujg {vc) statistic, we replace the ^^’s with 
(^i (x’c)’s and carry out the same computation. 

When estimating the covariance matrix from a finite number 
of mocks we introduce a bias which is corrected by multiplying the 
inverse covariance matrix by (1 — 79). 79 is defined by: 


79 = 


TXt, + 1 
ris - 1 


( 20 ) 


Here, rih is the number of bins used in the analysis and ris is the 
number of mocks used. 

To account for the effective volume of our sample, we scale 79 
by; 


r = 


2 Voverla 


VnGC + VsGC 


( 21 ) 


For the DRll sample, we have r = 0.49. 

The correlation matrix used when fitting the Ui (rc) statistic 
is presented in figure 


4 ANALYSIS 

4.1 Fiducial Cosmology 

Throughout our analysis we assume a fiducial ACDM cosmology 
with Qm = 0.274, Qb = 0.0457, h = 0.7, and Ug = 0.95. In this 
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bin 


Figure 1. The correlation matrix for uo^ {tc) computed for our analysis. 
The top figure depicts the covariance in ujq {tc) and the bottom figure de¬ 
picts the covariance in uj 2 (rc). Here, the Vc bins run from 46/i“^Mpc to 
198/i“^Mpc in bins of size 8/i“^Mpc. 


fiducial cosmology, the angular diameter distance to z = 0.57 is 
Da (0.57) = 1359.72 Mpc, the Hubble parameter is H (0.57) = 
93.56 kms“^Mpc“^, and the sound horizon is = 153.19 Mpc. 


4.2 Measuring the Correlation Function and Ui (tc) 


Our correlation functions are computed in radial bins of 8/i“^Mpc, 
running from 6/i“^Mpc to 198/i“^Mpc, and bins of 0.01 using 
the Landy-Szalay estimator ( |Landy & Szalay|1993| l: 


^ (r, li) 


DD (r, ii) — 2DR (r, fi) + RR (r, fi) 
RR (r, fi) 


( 22 ) 


and P (ko) = 20000/i“^ Mpc^ is the estimated power spectrum at 
the BAO scale. 

In addition to the FKP weight, each object is also assigned a 
systematic weight in order to account for various systematic errors 
in the data collection. The total systematic weight for an object, 
wtot, is defined in |Anderson et aL| ( [2014a| ) as: 

Wtot = {Wcp PWzf - 1) WstarWsee (24) 


Here, Wzf accounts for failures in acquiring the redshift of the near¬ 
est object and Wcp accounts for failures in acquiring redshift as a 
result of being in a close pair with a neighboring object. Further, 
Wstar and Wsee account for systematic relationships between the 
number density of observed galaxies, stellar densities, and seeing 
of data. 

Once we have computed and ^2 as described in the preced¬ 
ing section, we compute uq (tc) and 0 J 2 {tc) by applying Eq. 13 


using the filter defined in Eq. The integrals are computed over 
the fitting range and evaluated at the same re’s as the r bins in 


4.3 Fitting Methodology 

Our multipoles analysis closely follows the analysis conducted in 
[Anderson et al.| ( |2014b| ) and jXu et al.] ( |2012| ). We fit both the 
monopole and quadrupole simultaneously to our model for four 
nonlinear parameters: a, e, log and (3. This model is con¬ 

structed based on our expected observations given in Eqp^ 


^0,m (r) = ir) + ^0 (^) 


(25) 


6,m (r) = ^ 0 ,t (r) + A 2 (r) 


(26) 


where the subscript t denotes the templates given in Eq[^and the 
subscript m denotes our model. 

Our models for {vc) are derived from the templates for 
by applying Eq.pj]with the filter defined in Eq.I^to the expected 


observations for given in Eq. 18 In short, uji (rc) is a filtered 


integral of As in the multipoles fitting, when fitting uj^ (rc), we 
fit cjo {tc) and uj 2 (rc) simultaneously with the same nonlinear pa¬ 
rameters. 

We detail here the fitted parameters. Bq is a term that scales 
the monopole template. Before fitting, we estimate Bq from the 
offset between our model and the measured correlation function 
at r = 50/i“^Mpc. This estimation is then used to normalize the 
monopole to ensure that Bq ~ 1. To prevent unphysical negative 
values, we vary log (Bq) instead of Bq. We also apply a Gaussian 
prior centered at 0 with a standard deviation of 0.4 to keep Bq near 
1 . 

In order to improve the fit, we can choose to add polynomial 
terms Ao(r) and A 2 {r) to ujq (rc) and UJ 2 (rc) when fitting ujg (rc), 
similar to what is done when we fit the multipoles above. Ai (r) is 
defined by: 


where DD, DR, and RR are the number of galaxy-galaxy, galaxy- 
random, and random-random pairs separated by r and g. The ran¬ 
dom points are given by a set of randomly distributed points with 
the same selection function as our survey data. 

Each object is also given a EKP weight { [Eeldman, Kaiser &| 
|Peacock|1994j ) given by: 

1 3-n (zi) P (ko) 

where n (zi) is the number density at the redshift of the zth object 


A ( \ _ dig , dl^2 . 

Al (r) — —I—-—h ai ^3 (27) 

These polynomials are made up of linear nuisance terms used to 
marginalize out broadband effects like scale-dependent bias and 
our mismatch between the theoretical redshift-space distorted bi¬ 
ased galaxy correlation function and the true correlation function 
( jXu et al.|2012j ). These polynomial terms were used when working 
with the multipoles. However, they were not used with the uji (rc) 
statistic since we expect that polynomial terms are not needed when 
fitting uji (rc). 
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Tmin (/l~^Mpc) 

22 

30 

38 

46 

54 

62 

70 

78 

86 


1.008 ± 0.028 
0.998 ±0.014 
1.001 ±0.015 
1.003 ±0.015 
1.001 ±0.018 
1.000 ±0.017 
1.007 ±0.026 
0.991 ±0.013 
1.062 ± 0.030 


0.032 ± 0.022 
0.005 ± 0.015 
0.000 ±0.017 
0.001 ± 0.022 
0.000 ± 0.023 
0.000 ± 0.022 
-0.006 ± 0.040 
0.002 ± 0.028 
0.009 ± 0.083 


Table 1. Average and standard deviation for a and e for our 600 mocks as 
'f'min was varied. They were fit using (re) without polynomial terms. 


Our templates are based on the model for the correlation func¬ 
tion described in Sec. |3.3| with the streaming scale i n Eg . set 
to Eg = 1.4/i“^Mpc. The non-linear damping in Eq.[^is set to 
E^ = 6/i“^Mpc and E|| = ll/i“^Mpc. We also allow ^ to vary 
to better match the amplitude of the quadrupole. For /3, we apply 
a prior centered at 0.4 with a standard deviation of 0.2 to prevent 
unphysical values. Furthermore, since the e distribution is nearly 
Gaussian ( |Xu et al.|2012| ), we apply a 10% Gaussian prior on 1 ± e 
centered on zero to prevent unrealistic values for e. 

The non-linear parameters are fit using a simplex algorithm 
which steps through the parameter space until the • 



is minimized. Here, m is a column vector containing the model, 
dis 3. column vector containing the data and C is the covariance 
matrix defined in Sec. 13.41 

When we fit the multipoles, we adopt a fitting range of 
46/i“^Mpc to 200/i“^Mpc, which matches the range adopted in 


Anderson et al. (2014b i. For the {tc) fits, we analyzed the pos¬ 
sible fitting ranges by fitting our 600 mocks with varying Tmin^ 
and a fixed Vmax — 200/i“^Mpc for our fitting range and record¬ 
ing the best fit a and e values we obtained. Based on the results 
in Table[^ we concluded that an Tmin of 46/i“^Mpc would be the 
best for maintaining a close comparison between the two methods. 


4.4 Error Estimation 

Following the error analysis conducted in |Xu et a'L] ( |2012| ), the er¬ 
rors on both the uji (vc) and multipoles fits to data are estimated by 
computing the probability distribution, p (a, e), on a grid of (a, e). 
For each a and e in our grid, we fix a and e and fit the remainder 
of the parameters using as a goodness of fit indicator. Under the 
assumption that the likelihood is Gaussian: 

p (a, e) oc exp ^ ^ (29) 

and normalizing accordingly, we find that: 


p(a) = 

p (a, e) de 

(30) 

w.)= , 

f p (a, e) da 

(31) 


Since we assume Gaussian likelihoods, we can take the widths of 
the distribution a a and ae : 

p (a) {a — {a))‘^ da (32) 


Method 

a 

e 

(Xi (Xc) 

1.003 ±0.015 

0.001 ±0.017 

Multipoles 

1.004 ±0.016 

0.002 ±0.018 


Table 2. Average and standard deviation for a and e for our 600 mocks, 
fitted using cui (vc) without polynomial terms and the multipoles method. 


Method 

a 

e 

{Xc) 

1.029 ±0.020 

-0.003 ±0.023 

Multipoles 

1.025 ±0.014 

-0.010 ±0.019 

Anderson, et al. (2014b) 

1.025 ±0.014 

-0.010 ±0.019 


Table 3. Best fit values and standard deviation for a and e for DRll data, 
fitted using (vc) and the multipoles method. We also include as a refer¬ 
ence the fit results for multipoles with the de-wiggled template conducted 
in Anderson et al. (2014b). 


al = jp{e){e-{e)f de (33) 

as our error estimates. Here, {x) is the expected value of the distri¬ 
bution: 

For our error estimation, we range a from 0.7 to 1.3 inclusive 
with a spacing of 0.005 and e from -0.3 and 0.3 inclusive with a 
spacing of 0.01. This yields 121 a bins and 61 e bins. 


5 FITTING RESULTS 
5.1 Testing Mocks 

To appraise the performance of fitting the uji (vc) statistic, we ap¬ 
ply the two clustering estimators to the 600 PTHalos mocks and 
compare the results. Since the cosmology of the mocks is known, 
we apply it as the fiducial cosmology and use it in our fitting, we 
expect best fit values for a near 1 and best fit values for e near 0. 
The average and standard deviation for the fit parameters are shown 
in Table 1^ These results suggest that both methods give consistent 
values for a and e with discrepancies within our error bars. We 
note that the smaller errors returned by using the (xc) statistic 
are reflected in the x^ distribution from fitting the mocks. The x^ 
distribution produced by fitting (xc) yields a tighter distribution, 
centered closer to 1. 


5.2 Fitting Data 

We also test the performance of the two methods on data from 
DRll. Our results are summarized in TableHere, the errors 
on our parameters are computed using the analysis outlined in |4.4| 
Figure [^depicts our fit to the uj£ (xc) statistic for DRl 1 data. 

Our results for a are in agreement for the two methods tested 
here. The results we obtained for e also agree within 1 sigma. These 
numbers suggest that the two methods are consistent with each 
other. The results of our multipoles fitting are also in agreement 
with the numbers given in [Anderson et al.| ( |2014b| ), which is reas¬ 
suring. 


6 ADVANTAGES OF tj^rc) 


In this section we conduct a brief investigation into a possible ad¬ 
vantage of using the (xc) statistic. Based on Xu et al. (20101, we 
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(Mpc/h) 


Figure 2. Fits to DRll data using ui (vc). i = 0 is above and ^ = 2 is 
below. The crosses represent the DRll data and the blue line is our best fit. 

expect Ui (vc) to respond better to large scale abnormalities than 
the multipoles method. In order to test this response we conduct 
two tests comparing the performance of the U£ (tc) method against 
the multipoles method. 

6.1 Polynomial Addition of Extra Power 

Our first test consists of directly adding a polynomial to (^) to 
introduce abnormalities at large scales. The choice of polynomial is 
inspired in part by form of large scale power we observed in early 
BOSS DR9 data as well as in the quasar correlation function. The 
polynomial added, p (r), is of the following form: 



Here, x = rl27b and S' is a scale factor we adjust to change the 
amplitude. Figure [^depicts the addition of polynomials with vary¬ 
ing S' to ^ 0 - We apply the above treatment to the correlation func¬ 
tions for our 600 mocks and repeat the same fitting analysis applied 
earlier (Sec.|^. 

The best fit values for various scale factors S are presented 
in Table Again, since we know the fiducial cosmology is set to 
that of the mocks, we expect a to be approximately 1 and e to be 
approximately 0. 

We see that, as the value of S is increased, the average best 



Figure 3. An unmodified compared to 2 modified ^qs. 


Method 

a 

e 

Mult. S = 0 
Mult. S = 25 
Mult. S' = 45 
Mult. S = 65 

1.004 ±0.016 
1.004 ±0.016 
1.004 ±0.017 
1.004 ±0.017 

0.002 ± 0.018 
0.001 ±0.018 

-0.000 ± 0.020 

-0.001 ± 0.020 

u;£ (xc ), S = 0 
uj£ (xc ), S = 25 
lo£ (t’c) 5 S ^ 45 
uj£ {xc ), S = 65 

1.003 ±0.015 
1.003 ±0.015 
1.003 ±0.016 
1.003 ±0.016 

0.001 ±0.017 

0.000 ±0.018 
0.000 ±0.018 

-0.001 ±0.019 


Table 4. Average and standard deviation for a and e for our 600 mocks, fit 
using Ui (tc) (without polynomial terms) and the multipoles method with 
varying values of S. 

fit values and standard deviations for the multipoles fit remain rel¬ 
atively unchanged. The average best fit values and standard devia¬ 
tions for the Ui (xc) fit also remain relatively unchanged from the 
default of S' = 0. The U£{rc) fit performs similar to the multipoles 
fit in response to changes in S. 

Our results indicate that even without polynomial nuisance pa¬ 
rameters in the uj£ (xc) fit, it manages to match the performance of 
the multipoles fit [when extra large scale power is added], which 
uses polynomial nuisance parameters. 

6.2 Weight Application 

Our second test investigated how applying subsets of the weights 
(Wzf,Wcp,Wstar,Wsee) whcu Computing the correlation function 
for the DRl 1 data affected the performance of each fitting method. 
The SDSS-III BOSS analysis derived these weights in an effort 
to reduce the effect of observational systematics on data analysis. 
However, there are no errors computed for these weights. It would 
be advantageous if the ij£ (xc) method conducted without weights 
returned similar results to the uj£ (xc) method conducted with all 
the weights, making the absence of errors in the weight computa¬ 
tion less of a problem. The different weight applications, given by 
different forms of Eq.[^and their effect on ^£ (r) and uj£ (rc), are 
depicted in Figures]^ and 

Table|5]contains the results of our tests. We find that in the ma¬ 
jority of weight applications, although the multipoles fit performs 
slightly better, the two fitting methods remain comparable. As in 
the previous test, our results indicate that even without polynomial 
nuisance parameters in the uj£ (xc) fit, the cj£ (xc) fit nearly matches 
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Figure 4. The effects of different weight applications on computed from 
DRl 1 data. The top plot depicts and the bottom plot depicts ^2 • The top 
plot depicts the effects of different weight application on ^0 {f). 


Figure 5. The effects of different weight applications on (tc), computed 
from DRll data. The top plot depicts ojq (fc) and the bottom plot depicts 
a72 ( 7 * 0 ) • We see that the application of various systematic weights do not 
affect the shape of the (rc) statistic, making it a stable estimator in the 
face of systematics. 


Method and Weights 

Ao 

Ae 

Mult. All 

0.0000 

0.0000 

Mult. None 

0.0027 

0.0002 

Mult, wcp 

0.0043 

0.0016 

Mult. WZF 

0.0003 

0.0004 

Mult, {wqp + fuzf — 1) 

0.0011 

0.0006 

Mult. Wstar {wcp + WzF - 1) 

0.0004 

0.0003 

Mult. Wsee{wcp + WzF — 1) 

0.0006 

0.0009 

uji (rc); All 

0.0000 

0.0000 

(rc); None 

0.0037 

0.0024 

Wi (rc); WCP 

0.0066 

0.0012 

U£ (rc); WZF 

0.0004 

0.0019 

W£ (rc); {wcp + WZF — 1) 

0.0026 

0.0015 

LUi {rc); Wstar{wcp + WzF “ 1) 

0.0012 

0.0005 

Wi {rc); Wsee{wcp + WzF “ 1) 

0.0014 

0.0009 


Table 5. Absolute difference of best ht a and e for each weight application 
from the default of applying all weights. 


the performance of the multipoles fit in the face of systematics even 
though the multipoles fit uses polynomial nuisance parameters. 


7 CONCLUSIONS 

Observational systematics in cosmological surveys can impart sys¬ 
tematic fluctuations at large scales. It is useful to investigate new 
statistics that are less sensitive to these large scale fluctuations. In 
this paper, we consider one such statistic, Ui (rc). In order to de¬ 
termine the power of the uji (rc) statistic compared to the standard 
multipoles analysis, we applied both methods to fit for a and e in 
600 PTHaloes mocks and data from SDSS DRl 1. After comparing 
the best fit values and standard deviations returned by each method, 
we conclude that fitting with the Ui (rc) statistic without polyno¬ 
mial nuisance terms is comparable to the standard multipoles fit¬ 
ting. 

In addition, we also investigated the performance of each 
method in response to abnormalities in the correlation function at 
large scales. In line with the predictions in |Xu et aL] ( |2010| ) we find 
that fitting the uji (rc) statistic is less sensitive to abnormalities gen¬ 
erated by the addition of extra power (in the form of a polynomial) 
at large scales than the standard multipoles fit. The uji (rc) statis¬ 
tic matches the performance of fitting the multipoles even without 
the use of polynomial nuisance parameters. We observe similar re¬ 
sults when analyzing the response of the Ui (rc) statistic to changes 
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in the application of systematic weights to the data. Without using 
polynomial nuisance parameters, the (vc) statistic returns simi¬ 
lar results as those given by fitting the multipoles. 

The investigation here suggests that the Ui (vc) statistic could 
be advantageous compared to the standard multipoles method when 
dealing with abnormalities at large scales. Fitting the Ui (vc) statis¬ 
tic eliminates the need for the polynomial nuisance parameters 
which must be used in conventional methods for fitting the mul¬ 
tipoles. 
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